Diffusive redistribution of small spheres in crystallization of highly asymmetric binary 

hard-sphere mixtures 



O ■ 

<n : 

c : 
s . 

o : 



o 



O 

o 



> 
o 
in 

On 
O 



X 



Wen-Sheng Xu, Zhao-Yan Sun*, and Li-Jia An* 
State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, 
Chinese Academy of Sciences, Changchun 130022, People's Republic of China 

(Dated: June 13, 2011) 

We report a molecular dynamics study of crystallization in highly asymmetric binary hard-sphere 
mixtures, in which the large spheres can form a crystal phase while the small ones remain disordered 
during the crystallization process of the large spheres. By taking advantage of assisting crystal nu- 
cleation with a patterned substrate, direct evidence is presented that there is a close link between 
the diffusive redistribution of the small spheres and the crystal formation of the large spheres. 
Although the addition of a second component with large size disparity will not alter the crystal 
structure formed by the large spheres, the density profile of the small spheres displays correspond- 
ing changes at different crystallization stages and closely relates to the crystal growth, suggesting 
possible effect of small spheres on the crystallization kinetics. 



Crystallization of hard-sphere colloids has been the 
subject of intense research since assemblies of hard 
spheres constitute a very simple model but have poten- 
tial important applications and exhibit complex behav- 
ior. Pioneering work by Alder and Wainwright in 1957 [1] 
demonstrated that the monodisperse hard-sphere sys- 
tem crystallizes at high enough volume fraction, driven 
by purely entropic effects. Over the last few decades, 
understanding for crystallization of hard-sphere colloids 
has been greatly improved due to the development of 
new experimental techniques and novel numerical 

methods [5, 6]. The binary hard-sphere mixtures with 
moderate size disparity also attract considerable interest 
for both practical and fundamental reasons, e.g., vari- 
ous types of binary crystals can be formed in such sys- 
tem [7H10| and a nonstandard nucleation and growth sce- 
nario has been proposed by studying binary hard-sphere 
mixtures with a small to large sphere diameter ratio of 
q = u s /(Ti = 0.9 (<r s and ai denote small and large sphere 
diameters) at a total volume fraction of 0.58 In 
contrast, much less attention was received on the crys- 
tallization of highly asymmetric binary hard-sphere mix- 
tures, in which the large spheres can form a crystal phase 
while the small ones remain disordered. The main rea- 
son resides in the fact that the addition of a much smaller 
second component will intuitively not affect the crystal 
structure and hence properties of the resulting material. 
However, the crystallization kinetics in such system is 
still unclear. 

In this work, we focus on the crystallization of binary 
hard spheres with large size disparity. In such system, 
the presence of the small spheres can induce attractive 
forces between the large spheres, called the depletion ef- 
fect 12[ , which is most clearly present in colloid-polymer 
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mixtures. By taking this effect into account (through ei- 
ther the effective one-component model treatment or a di- 
rect study of true binary mixtures) , extremely rich phase 
behavior, including a stable isostructural solid-solid tran- 
sition, a metastable fluid-fluid transition and two types 
of glassy states, has been revealed in the highly asym- 
metric binary hard-sphere system (l3| . However, less is 
known about the role played by the smaller component in 
nonequilibrium processes such as crystallization process. 
In particular, is there a possible link between some prop- 
erties of the small spheres and the crystal formation of 
the large spheres? In a recent paper on crystallization of 
colloid-polymer mixtures 14[, anomalous crystallization 
kinetics was observed, which was postulated to be due to 
the diffusive redistribution of the disordered polymers. 
Our work can also test this postulation since the highly 
asymmetric binary hard-sphere mixtures can serve as a 
reference model for mixtures of colloids and polymers. 

In principle, one can study these questions from the 
homogeneous crystallization of binary mixtures. How- 
ever, because the nuclei are formed randomly throughout 
the system, it is impossible to control the orientation of 
the growing crystal and the shape of the forming crystal- 
lite is often irregular. This makes it difficult to charac- 
ter the properties of the small spheres during the crystal 
formation. To alleviate this, a patterned substrate can 
be introduced into the system as the seed of crystalliza- 
tion, then the direction of the growing crystal can be 
controlled by the substrate and the crystal grows almost 
layer-by-layer. Although this also introduces a degree 
of complexity compared to the homogeneous crystalliza- 
tion, we expect that the conclusion in this work can be 
generalized to the case of homogeneous crystallization. 
We find that the addition of a much smaller second com- 
ponent will not alter the crystal structure formed by the 
large spheres. However, by monitoring the density pro- 
files of both species during the crystallization of the large 
spheres, we observe that the density profile for the small 
spheres displays corresponding changes at different crys- 
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tallization stages and closely relates to the crystal growth 
of the large spheres. 

We use a (100) patterned substrate and the lattice 
spacing is a = 1.576(7;. The patterned substrate is set in 
the x-y plane and located in the middle of z-direction of 
the simulation box. Periodic boundary conditions are ap- 
plied in all three directions. The spheres in the substrate 
are fixed, i.e., when the other spheres collide with those 
in the substrate, they bound off as if the spheres in the 
substrate have infinite mass. The previous work 15, 1(| 
has shown that the forming crystal for the monodisperse 
hard spheres has a purely face-centered cubic (fee) struc- 
ture under the influence of this patterned substrate. For 
binary hard spheres with large size disparity, the struc- 
ture of crystals formed by large spheres is still fee, which 
will be discussed later. All binary hard-sphere liquids 
are quickly compressed to a large sphere volume frac- 
tion of = irNiaf / &V = 0.535 (well within the coex- 
istence region of hard spheres, AT; stands for the large 
sphere number and V the volume of the simulation box) 
by using the Lubachevsky-Stillinger (LS) algorithm [17 1 
with a compression rate of 0.02. Crystallization process 
is then simulated using a collision-driven molecular dy- 
namics simulation technique In this work, we only 
focus on the crystallization of binary mixtures with con- 
siderably large size ratios, and therefore the size ratio q 
ranges from 0.1 to 0.2. The volume fraction of the small 
spheres cf> s — irN s a^ /6V with N s the small sphere num- 
ber, covers the range from 0.001 to 0.009. The large and 
small spheres have the same mass m. Since supersatu- 
ration depends on the composition of binary mixtures, 
usage of the same compression rate will result in differ- 
ent supercoolings of the initial liquids. This may affect 
the speed of crystallization and further discussion will be 
given later. The temperature T is set to be 1. Length 
and time are given in units of 07 and ai \fmjkj^T with ks 
the Boltzmann constant. The number of large spheres is 
fixed at Ni = N bulk + N sub = 1000 with N sub = 100 and 
Nbuik = 900, where N sub and N bu ik are the sphere num- 
bers in the substrate and in the bulk. The corresponding 
dimensions of the simulation box are L x = L y = 7.88(7; 
and L z = 15.76(7/ . Although AT; = 1000 may seem too 
small a number to perform simulations of hard-sphere 
crystallization, we note that (i) the maximum value of 
N s is 16822 due to the small size ratio, (ii) our aim 
is to identify the possible role of small spheres during 
crystallization process, rather than to focus on the form- 
ing crystal by the large spheres, and (iii) results on the 
large spheres are qualitatively consistent with the previ- 
ous work [15] ]. in which a much larger system was used 
for the same patterned substrate. 

We use a spherical harmonic method 19|, 2(| to identify 
the large spheres in a crystalline environment. It defines 
particles as being ordered in terms of coherence of the 
spherical harmonics of a particle with its neighbors. In 
our work, the neighbors of the large sphere are defined 
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FIG. 1: Sample-averaged distribution of the local bond- 
orientation order parameter of the large spheres at varying 
time for q = 0.1 and (f> s = 0.005. The peak emerging at 
around 0.57 indicates the formation of a purely fee crystal 
structure. The result is averaged over 200 independent runs. 



as those within a cutoff distance of 1.4(7; from the target 
sphere. We use the criterions with d c — 0.5 (This cri- 
terion identifies candidate crystalline particles as having 
d 6 > d c , please see Ref. [H, |2(| for the definition of de) 
and N c = 10 (This criterion requires the target particle 
to have at least N c neighbors that also satisfy d§ > d c ). 

We first present the distribution of the local bond- 
orientation order parameter for the large spheres at 
q = 0.1 and <f> s = 0.005 in Fig. 1 (see Supplemen- 
tary Video SI for a full time series). The local bond- 
orientation order parameter is defined as Qe.iocai — 
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spherical harmonics with I 
number of neighbors for the jth sphere. The values of 
Qe.iocai for perfect structures of fee, hexagonal closed 
packed (hep), and the body-centered cubic (bec) crys- 
tals are 0.575, 0.485, 0.511, respectively. It is clearly 
seen in Fig. 1 that a peak at around 0.57, which indi- 
cates the formation of the fee structure, emerges and its 
height increases as crystallization proceeds. No apparent 
peaks appear at other positions, thus the forming crys- 
tal structure under the influence of (100) patterned sub- 
strate ispurely fee, which is consistent with the previous 
work |l5l . [l6| . This indicates that the crystal structure 
formed by the large spheres is not altered by the addition 
of very small spheres. 

We show the diffusion behavior of both species in Fig. 
2 by calculating the mean squared displacements (MSD). 
We observe that the ballistic motion (< Ar 2 (i) >^ t 2 ) 
at short times crosses into the diffusive behavior (< 
Ar 2 (t) >~ t) at long times for the small spheres. For 
the large spheres, the diffusive region cannot be reached 
since the diffusion ceases when the crystallization com- 
pletes (the crystal fraction reaches 1). At that time, a 
typical large sphere has diffused only about 1 large sphere 
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FIG. 2: Sample-averaged mean squared displacements (MSD) 
of both species at (a) varying 4> s with q — 0.1 and (b) vary- 
ing q with (j> s = 0.005. N a denotes small sphere number or 
large sphere number in the bulk. The red dash lines with 
slopes 1 and 2 have been included to stress the fact that the 
small spheres undergo typical diffusion of simple liquids. The 
results are averaged over 50 independent runs. Insets: ex- 
tracted slopes of MSD for both species. 



diameter. This is about 10 times smaller than that in the 
homogeneous crystallization plj . indicating the high ef- 
ficiency of the substrate in directing crystal growth. To 
better character the diffusion behavior of both species at 
different q and cj> s , we also present the extracted long- 
time slopes of MSD in the insets. For the large spheres, 
the slopes are obtained by extracting long-time data be- 
fore the cessation of diffusion. At fixed q, the slopes for 
both species decrease as <p s increases (Fig. 2(a)). Surpris- 
ingly, the slope for large spheres increases as q increases 
at fixed <p s (Fig. 2(b)), which indicates that the speed of 
crystallization is fast for large size ratio at fixed volume 
fraction of small spheres. This is also confirmed by the 
growth speed of crystal fraction (data not shown). We 
speculate that this is because of a change in the driv- 
ing force since the supersaturation of the initial super- 
cooled binary liquids is different for various mixtures in 
our study. As the small spheres undergoes diffusion of 
typical simple liquids, one may think that the density 
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FIG. 3: Sample-averaged density profiles of the large (upper 
panel) and small (lower panel) spheres for q — 0.1 and (f> s — 
0.005 at (a) t = 100 and (b) t = 1000. The yellow dash lines 
indicate the substrate position. H(t) denotes the thickness of 
the crystalline film (see text for its definition) and z a (t) the 
position relative to the substrate where the maximum local 
density for small spheres exhibits. The density profiles are 
averaged over 200 independent runs. 
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FIG. 4: Connection between the position where the maximum 
local density for the small spheres displays and the thickness 
of the crystal film formed by the large spheres at q = 0.1 and 
d>„ = 0.005. 



profile of small spheres during crystallization process is 
fully random. However, we will show below that the dis- 
tribution of the small spheres has a close connection with 
the crystal formation of the large spheres. 

In Fig. 3, we present the density profiles of both species 
at t = 100 (Fig. 3(a)) and t = 1000 (Fig. 3(b)) (see Sup- 
plementary Video SI for a full time series). The density 
profiles, denned as the average number of spheres per 
unit volume along the z-axis, are accumulated in bins 
of width Sz = O.lcr;. Obviously, there exists a correla- 
tion between the density profiles for both species. For 
instance, at the stage of crystal growth (t — 100), a max- 
imum local density for the small spheres is found at the 
position where the thickness of the crystal formed by the 
large spheres can be roughly identified (The detailed dis- 
cussion will be shown later). While, the small spheres will 
distribute uniformly in the voids of the resulting crystal 
after crystallization completes (t — 1000). 

The more clear support for connections between the 
small sphere distribution and the large sphere crystal 
formation is provided by calculating the maximum lo- 
cal density of the small spheres pT ax (t) an d the thick- 
ness of the crystalline film H{t) for the large spheres at 
time t. p™ ax {t) is determined by searching for the max- 
imum local density of the small spheres, which can be 
readily identified before the complete of crystallization, 
and the corresponding position relative to the substrate 
is denoted as z s (t). We determine H(t) as the last peak 
position relative to the substrate in the large sphere den- 
sity profile, where the peaks are defined as local densities 
larger than pi+p t , where pi indicates the average density 
of the large spheres along the z-axis and pt is a thresh- 
old. We use p t — 1.0 and the qualitative result is not 
dependent on the choice of p t . The resulting H(t) and 
z s (t) are shown in Fig. 4. As can be seen, they have 
similar behavior with t, i.e., the growth of the crystalline 
film is accompanied by a shift of z s (t) to large distances 
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relative to the substrate before crystallization completes. 
H(t) saturates at 7.05o7 and z s (t) varies randomly when 
crystallization completes. Thus, the results unambigu- 
ously demonstrate a coupling between the distribution 
of the small spheres and the crystal growth of the large 
spheres during the crystallization process. In Fig. 5, 
we also show the maximum local density for the small 
spheres p™ ax and its connection to the growth of the 
large sphere crystal fraction during crystallization pro- 
cess. We observe that (i) p n s lax has a large value in the 
initial supercooled liquid and the corresponding position 
is near the substrate (see Fig. 4), which is due to the pres- 
ence of the substrate during the compression; (ii) when 
entering into the stage of crystal growth (t > 1), p" lax 
at first drops from ~ 12.2 to ~ 11.4 and then remains 
nearly constant at ~ 11.4 for t < 100. The latter ob- 
servation suggests that there is first an expulsion of the 
small spheres followed by a diffusion back-in during the 
crystallization of the large spheres. In other words, the 
small spheres diffusively redistribute during the crystal- 
lization process; (iii) on the approach to the complete 
of crystallization, p™ ax dramatically increases to ~ 12.0 
followed by a fast decrease, indicating the non-balance 
between the expulsion and back-in of the small spheres 
at this stage; (iv) after the complete of crystallization, 
the small spheres uniformly distribute in the voids of the 
crystal formed by the large spheres, thus p™ ax remains 
unchanged at about 10.5. The result again clearly shows 
that the distribution of the small spheres exhibit corre- 
sponding changes at different crystallization stages. And 
the observation for the small spheres may also help to 
better understand and distinguish the different crystal- 
lization stages. 

In summary, by taking advantage of assisting crystal 
nucleation with a patterned substrate, we have demon- 
strated very directly a close link between the density 
profile of the small spheres and the crystallization of 
the large spheres in the highly asymmetric binary hard- 
sphere mixtures. Our finding may potentially explain the 
experimentally observed anomalous crystallization kinet- 
ics in colloid-polymer mixtures. However, since our work 
only gives evidence that the small spheres do play some 
role in the crystallization process of the large spheres 
and the crystallization kinetics (such as the growth rate 
of the crystallite size) cannot be probed due to the small 
size in our simulations, we hope that our finding will 
trigger relevant experiments on highly asymmetric bi- 
nary hard spheres, in which crystallization kinetics can 
be measured. 
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FIG. 5: Connection between the maximum local density for 
the small spheres (upper panel) and the growth speed of crys- 
tal fraction for the large spheres (lower panel) at q — 0.1 and 
cj>s = 0.005. 
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